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ABSTRACT 

The helioseimic investigation shows that the convective overshooting can penetrate 
0.37Hp to the location where the temperature is 2.5 x 10 6 K which is the typical 
temperature of the reaction Li 7 (p, a)He A . This indicates that the overshooting mixing 
should be involved in investigating the solar Li abundance problem. Observations of 
Li abundance of solar-twins show that the sun is not very peculiar. The overshooting 
mixing should be also involved in investigating Li abundance in clusters. However, the 
fully overshooting mixing with the length of 0.37f/p results in too much Li depletion to 
fit the observation in the solar case. Therefore, using the diffusive process to describe 
the overshooting is more suitable. The diffusive overshooting approach requires the 
turbulent r.m.s. velocity in the overshooting region to calculate the diffusion coefficient. 
Turbulent convection models(TCMs), which are suggested by helioseimic investigation, 
can provide the turbulent properties in the overshooting region. However, TCMs are 
often too complex to be applied in the calculations of the stellar evolution. It is an 
easier way to use the asymptotic solution of TCMs. In this paper, I use the asymptotic 
solution of Li & Yang's TCM, which results in agreements in both solar sound speed 
and solar Li abundance, to investigate Li abundance in clusters (Hyades, Praesepe, 
NGC6633, NGC752, NGC3680 & M67). It is found that the overshooting mixing 
leads to significant Li depletion in the old clusters^ > lGyr) and little effect in the 
young clusters(i < lGyr). 

Key words: convection - diffusion - stars: abundance - galaxies: clusters: individual: 
Hyades, Praesepe, NGC6633, NGC752, NGC3680 & M67. 



1 INTRODUCTION 

Li is a fragile element which can be burned via the reaction 
Li 7 (p,a)He 4 at the typical temperature T w 2.5 x 10 6 K. 
Therefore, the abundance of Li at the stellar surface is a 
good tool to probe the property of the mixing in the stellar 
interior. Li abundance can be derived by the spectroscopic 
analysis. Abundant observations of Li abundance are accu- 
mulated in the recent few decades. 

The observations of open clusters show some gen- 
eral relations betwe en Li abundance and T e //(see e.g., 
IXiong fc Dend (2009) and references therein). There are 
three remarkable properties: i) Li abundance decreases when 
T e j f decreases, ii) there is a Li gap near T e f / ~ 6500AT in the 
clusters with the age no younger than the Hyades's, iii) the 
observations show dispersions of Li abundance which can't 
be thought as the observational uncertainties. The second 
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could be due to th e atom i c diffusion and radiative a c celera - 
tion(i.e..lMichaudl Jl986l) : iMichaud fc Charbonneaul l|l99ll) ; 
IXiong fc Dend (|2009l )L 



There are some inconsistencies between the obser- 
vations and the relation of Li abundance vs. T e ff pre- 
dicted by the standard stellar evolution model. With the 
mixing length parameter of solar calibr ation and the old 

solar c hemical composition (Z w 0.02, IGrevesse" fc Noels! 

<|l993l ); IGrevesse fc Sauvall (|l998l )) which results in cor- 
rect location of the base of the solar convection zone, 
the standard stellar evolution model leads to too much 
Li depletion during the pre-main sequence( PMS) stage to 
fit the observations of young open clusters ( | Ventura et al.l 
Il998bl : iPiau fc Turck-Chiezd 12002ft . On the other hand, 
based on the standard stellar evolution model, which 
uses the convection as the only process of mixing, it 
has been found that Li depletion almost doesn't occur 
during main sequence(MS) stage for solar mass stars. 
This fails to explain the observations especially for the 
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old clusters in the low-temperature rang, and there 



MSfD'Antona & Mazzitclli 1984: Proffitt & Michaud 


198S 


D'Antona & Mazzitellil 199 


4: Piau & Turck-Chiezc 


2002 


D'Antona & Montalban 2003). Richard et al.l (2002 


2005 



have investigated low-metallicity stars and have found that 
extra turbulent mixing is required to offset the atomic dif- 
fusion. 

There are some possible mechanisms to de- 
plete Li in MS: mass- lossdHobbs et alj 1 19891 : 
ISchramm e t al. 1990; Swcnso n fc, Faulknerll 19921 ). rotational 
mixin g jpinsonneault et alj 1990l . 



1999: Charbormel et alj 



19921. 1 1994). internal waves mix ing |Montalbanl 1 1994 



Montalban 



Schatzmaiilll996l . [2000) 



lent mixing (Michaud 1986; 
Il99ll : iRichard et alT 120021. 



diffusion and turbu- 
Michaud fc Charbonneaul 
120051), and overshoot- 



ing mixingdStraus et alj Il976l: ISchlattl fc Weiss! 1 19991 : 
IXiong fc Denell2009l ). ISestito fc Randichl (|2005l ) have stud 



ied the evolution of Li abundance and have concluded 
that "gravity waves can be excluded as the main agent 
responsible for MS Li depletion, while slow mixing induced 
by rotation might explain to some extent the empirical 
behavior" . The rotational mixing, which is the most popular 
extra-mixing model, might reproduce the observational 
properties of the cluster Li abundanc e not only on the 
deplet ion but also on the dispersionfe.g..|Pinsonneault et all 
(1999)). Recently. IChristensen-Dalsgaard et alj ( 2011 ) have 
found that the convective overshooting could penetrate 
0.37Hp(Hp is the local pressure scale height) to the 
location where T = 2.5 x lO 6 -^. This indicates that the 
overshooting mixing should be involved in studying the solar 
Li problem. Investigations on Li abundance of solar-twins 
shows that the solar Li abundance seems to be not very 



differe n t comparing with sol a r- twin s ( see e.g., King et all 



ll 19971): iMelendez fc Ramfrej ||2007D : iLubin et al. I l|20ld) : 



ICastro et alj |201ll )). Therefore, the overshooting mixing 
should be also involved in investigating the Li abundance 
in open clu sters. Using Xio n g's(19 85) theory of non-local 
convection, IXiong fc Dend |2009l ) have found that Li 
is significantly depleted during the MS stage when the 
overshooting mixing is present. 

Recent helioseismic investigations have suggested 
using turbulent co nvection models(TCMs) to dea l wit h 
the overshooting (|Christensen-Dalsgaard et all 1201 if ). 
The TC Ms are based on fully hydrodynamic moment 
equa t ionsllXiongl Il985l: IXiong fc Den3l200ll: ICanutd Il997l 
19981 , Il999l : iDeng et alj|2006l : iLi fc Yandl2007l : lYang fcU 
20071 ) . Therefore, they are more physically reasonable 
than the phe nomenologic al theories in the framework of 
mixing length (|Hinzdll975l 'l. However, the TCMs are highly 
non-linear and are hard to be applied in stellar struc- 
ture and evoluti o n Based on t he TC Mj Li fc Yangj 120071 : 
lYang fc Lil 120071 ). I Zhang fc Lil (|2012bT ) have investigated 
the overshooting region and have found the asymptotical 
solution of the turbulent fluctuations in the overshooting 
region. Using the diff usive overshooting a pproach with the 
asymptotical solution, IZhang fc Lil l|2012al ) have found that 
the solar sound speed and the solar Li abundance can be in 
agreement with the observations. 

In this paper, I use the diffusive overshooting approach 
and the asymptotical solution of the Li & Yang's(2007) 
TCM to investigate the Li abundance in six open clusters: 



Hyades, Praesepe, NGC6633, NGC752, NGC3680 & M67. 
The stellar modeling method is described in Section 2. The 
diffusive overshooting approach is introduced in Section 3. 
The numerical results are shown in Section 4. Main conclu- 
sions and discussions are in Section 5. 



2 STELLAR MODELING 

In order to investigating Li abundance in clusters, I calculate 
the evolution series of the stellar models with different mass, 
metallicity and age. Six clusters are taken into account: three 
intermediate ones(Hyades, Praesepe & NGC6633) and three 
old ones(NGC752, M67 & NGC3680). 

The modified stellar evolution code (|Paczvnskil Il969l ). 
which has been originally described by Paczynski and Ko- 
zlowski and has been updated by Sienkiewicz, is used to 
calculate the ste llar evolution mode ls. The OPAL equation 
of state tables jRogers et all Il996l) are used. The OPAL 
opacity tables (|lglesias fc Rogers! 1 19961 ) are used in the re- 
gion of IgT > 3.95, and t he Alexander's opacity tables 



^Alexander fc Fergusonll994 ) are used in the r egion of IgT < 
3.95. The nuclear reaction rates are from iBahcall et al.l 
jl995l ) in the calculations. 

The metal mixture of all clusters is assumed as the solar 
metal mixture (|Grevesse fc Sauvallll998l ). The metallicity of 
each cluster is calculated as: 



[Fe/H] 



z = 10 



where (-fO© 



<#> * 



(i) 



0.023 is the ratio of th e metallicity to the 
hydr ogen abundance at the solar surface (|Grevesse fc Sauvall 
1998), [Fe/H] the metallicity dex of corresponding cluster, 
and X = 0.7 the initial hydrogen abundance of all stel- 
lar models. The solar comp osition used in t his pa per is not 
the latest observation(e.g., lAsplund et all (120091 ) ), for the 
reason that the new observational composition results in in- 
consistency on the location of the base of the solar con- 
vecti on zone between the solar model and the helioseismic 
datadBasu fc Antial l2004: IBahcall et al.l [20051: lYang fc Bil 



12007! : IChristensen-Dalsgaard et al.1 120091 : iBi et al j l201ll ) 

However, the location of the base of the convective en- 
velope(BCE) is crucial for calculating the Li abundance, 
since Li burning speed is mainly determined by the tem- 
perature of the BCE. T herefore, the solar composition of 
iGrevesse fc Sauvall {l998), which leads to the correct convec - 
tive boundary in th e solar modelfe.g.. IBahcall et all {2005): 
iPaxton et all (|201ll )). is adopted. 

The metallicity dex, age and the sources of Li abun- 
dance observations of concerned clusters are listed in Table 
1. The mass ranges of stellar models of each cluster are: 
0.80 ~ 1.20M© for Hyades, Praesepe, NGC752 and M67, 
0.70 ~ 1.2OM for NGC6633 and NGC3680. The mass-step 
is AM = O.O5M . 

The convection zone is assumed to be fully mixed. The 
convective boundary is defined by the Schwarzschild criteria. 
The convective heat flux is calculated by using the mixing 
length theory with the parameter a — 2.08 according to the 
calibration of the solar model by using those input physics. 

All of the stellar model series evolve from the PMS stage 
with the center temperature Tc = lO 6 -^ to the age of cor- 
responding cluster(Table 1). The time step is no more than 
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Table 1. The age and metallicity dex of clusters. The ref- 
erences of the age and Li abun d ance o f eac h cluster are as 
follows : fl lHobbs fe Pilachowskl jl986ft: (2lSoderblom et alj 
||1993eJ|; OlSoderblom et al.l jl993bh; (4lBalachandranl <1995tl; 
(5lDinescu et al.l jl995h; ( 6Uefiriesl Jl997f); (7lPerrvman et al l 



jl998h; (81 Jones et aD ll 199911 ; OlRandich et al.1 
' Pasquini et al l feOOl); (111 Jeffries et al.1 
Sestito et alj (20o3 ): (13 fearbara et al.l l|2009l V 
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Figure 1. The evolutional history of surface Li abundance, center 
temperature Tq and the BCE temperature TgoE of the solar 
model with fully overshooting mixing in 0.37-ffp. 



0.5 percent of the age of corresponding cluster. All of the 
stellar models include more than 2000 mesh points. 



3 THE CALCULATION OF LI ABUNDANCE 
AND THE OVERSHOOTING APPROACH 

3.1 The diffusive overshooting mixing approach 

The equation of the evolution of Li abundance in the stellar 
interior is as follows: 
BY 
~dt 



[(4-7rr pj D- 



dM r 



dM r 



pRXY 



(2) 



where Y is Li abundance fraction, X the hydrogen abun- 
dance fractio n, R the ratio of the reac tion Li 7 (p,a)He 4 
calculated by ICaughlan fc Fowled (|l988T l and D the diffu- 
sion coefficient due to the convection and the overshoot- 
ing. The convection zone is assumed as fully mixed, thus 
D = +00, (Vr ^ Vad), where V_R,V a d are the radiative 
temperature gradient and the adiabatic one. 

Li abundance is calculated in two cases: i) no overshoot- 
ing, ii) taking the downward overshooting of the BCE into 
account. If the overshooting is absent, D = 0, (Vr < V a d). If 
the overshooting is present, D — Dov, (Vij < V Q d), where 
Dov is the diffusion coefficient of the overshooting mixing. 

If the overshooting mixing is set to be instantaneous 
mixed as widely adopted in traditional, the sun shows ex- 
cessive Li depletion. A fully mixed overshooting region w ith 
the length of 0.37ffp (|Christensen-Dalsgaard et al.ll201ll ) re- 
sults in A(Li) « —4 as sho w in Fig.l. It is too l ow to fit the 
observation A(Li)@ = 1.05 (|Asplund et alJl2009T ). 

The fully overshooting mixing is excessively strong. The 
weak overshooting mixing described as a diffusive process 
is more suitable. The diffusion c oefficient of oversho oting 
mixing Dov is assumed as follows ! Zhang fc Lill2012al ): 



Dov = CxHpVk 



(3) 



where Cx i s a par ameter and Cx ~ 10 a ccording to 
I Zhang fc Li (|2012ah that is in consistent with iDeng et al.1 
l|l996l )(Eqs.(27') & (29)), k the turbulent kinetic energy in 



the overshooting region. 

The turbulent kinetic energy k can 
by solving TCMs which are suggested 



be obtained 
by the he- 



lioseismic investigation l|Christensen-Dalsgaard et al.1 l201ll ) . 
TCMs sho w that fc exponent i ally deceases in the overshoot- 
ing regiondXiond 11989 . 11989; IDeng et al.1 120061 : IZhang fc Lil 
2i Mi!). |2012bl ). This is in agr e ement with the numerical sim- 
ulations^. g., iFrevtag et al.1 |l996)). It is a good choice to 
solve TCMs numerically to get the profile of fc in the over- 
shooting region, but that is too complex. An easy way is 
to use the asymptotic solutio n of TCMs in the ov ershooting 
region. In the previous work dZhang fc Ll 2012bl). we have 
investigated the TCM developed by iLi fc Yand l|2007h and 
have obtained the asymptotic solution in overshooting re- 
gion as follows: 

k^k c (^-) e ,(Pe>i) (4) 



-Pc' 



and 



fc = 0, (Pe < 1) 

In Eqs.(4) and (5), Pe = pCpivT 



(5) 



— — is Peclet number, Pc the 
pressure at the BCE, 9 the exponential decreasing index of 
the turbulent kinetic energy fc, and fcc the turbulent kinetic 
energy at the BCE. 

It is crucial to get the turbulent kinetic energy at con- 
vective boundary (i.e., fcc) which appears in Eq.(4). A local 
convection theory(e.g., the mixing length theory) leads to 
fc = at the convection boundary. In the TCM, the non- 
local effect (turbulent diffusion of fc) is taken into account. 
In the framework of TCM of a convective zone it is possi- 
ble to define a point 'B' separating two regions shown in 
Fig. 2: the first one where the kinetic energy profile is well 
described by the local TCM: fcjvi ~ fcL (with fcjvz, and fcL 
the non-local and local kinetic energy respectively); and the 
second one between the point 'B' and the boundary of the 
convective region 'C where the non local terms dominate 
and d In fcjvL /d In P = 6 (with 8 a constant determined only 
by the parameters of the TCM). As a consequence, the log- 
arithmic value of the kinetic energy at the boundary of the 
convective zone (Inkc) can be obtained by linear extrapola- 
tion from the value of InkL at point 'B' with a slope 8. The 
profile of InkNL in the convection zone is the combination 
of the dashed line and the left part of the solid line to point 
B. Point 'B' is the the location where the turbulent diffu- 
sion becomes significant, so that dlnfci/dlnP = 9 at 'B'. 
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Convection zone 



The deep convection zone 



Convective 
boundary (C) 




-Ink, 



B: the location of dlnk/dlnP=e 



In P 



Figure 2. A sketch helps to understand the estimate of kc in 
non-local convection theory. Point 'B' is the location where the 
turbulent diffusion becomes significant, k^ is turbulent kinetic 
energy determined by the localized TCM, and /cjvl is turbulent 
kinetic energy determined by the non-local TCM. 



This proper ty results in t he max imum of kc(the maximum 
of diffusion, I Zhang fe Lil (|2012bl ) ) . 

According to the maximum of di ffusion, kc can be es- 
timated as follows I Zhang fc Li|[2012bl '): 



3 11 

k c = —k^ 



1 ai/3g\{V R - V ad ) 
e pc P 



(6) 



where j3 = — (dlnp/dlnT)p is the expansion coefficient, g 
the gravity acceleration, A = (4acT 3 )/(3«;p) the thermal 
conduction coefficien t and the location of point 'B' can be 
calculated as follows (|Zhang fc Li|[2012bl '): 



4C s o;c I 



= 1 



(7) 



where tb and rc are the radius of B and the convective 
boundary respectively, / = ot\Hp, ai,C s , Ck are parameters 
of the TCM and uj c = l/(3C fc ) + 1/3. 

The para meters of the TCM can be derived by ex- 
periments( e.g., Gibson & Launder (1976)), analyzing TCM 
itself (e.g., ICanuto fc Dubovikovl (Il998h) and the calibra- 
tion of the solar model (|Zhang fc Lill2012aT ). The appropri- 
ate paramet e r set o f the TCM for the solar case is listed in 
IZhang fc~L i (2012a). Accordingly, the value of the param- 
eters involved in this paper are: 6 — —4, a?i = 0.93, C s = 
0.08, C k = 2.5. 

Equation (4) is similar to Eq.(4) in Ventura et 
al.'s(1998) paper, since both of them describe the kinetic 
energy as a power law of the pressure in the overshooting re- 
gion. The difference between Ventura et al.'s(1998) descrip- 
tion and this paper is the calculation of the initial turbulent 
velocity. In Ventura et al.'s(1998) paper, it is determined by 
extrapolating the turbulent velocity distribution function re- 
sulting from the local convection theory near the convective 
boundary. In this paper, it is calculated based on the TCM 
which is a non-local convection theory. 



3.2 The initial Li abundance at ZAMS 

In the PMS stage, the stellar activity is significant thus 
the magnetic field could affect the stellar structure and 
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Figure 3. Li abundance of the young cluster Pleiades. Observa- 
tion data: circles and triangles are the observational Li abundance 
and the upper limits respectively. The solid line is the profile of 
Eq.(8), and presents the approximate profile of the low envelop of 
the Pleiades Li abundance. Dashed lines are the ZAMS Li abun- 
dance predicted by standard stellar model in which the PMS Li 
depletion is taken into account. From left to right, the dashed lines 
correspond to the Hyades, Praesepe, Pleiades case, respectively. 



the depletion of Li. Abundant observations show that 
the standard stellar theory underestimate the radii and 
overestimate the effec ti ve temperature of l ow mass stars 
dTorres fc Ribasl 120021: iBerger et all 120061: iMorales et all 
l200Sl , l2010r ). lTorres et all l|2006h and lRibaj l|2006l ) have sug- 
gested that those inconsistencies may result from the stellar 
activity because the surface activity and the magnetic field 
reduce the efficiency of convection and then the star should 
increase the radius to produce the required total radiative 
flux. With the mixing length parameter resulting from the 
solar calibration and the o ld solar chemical composit ion(e.g., 
iGrevesse fc Noels! (ll993r);lGrevesse fc Sauvall lll998l)'). it has 



been found that! Ventura et all 1998bl : iPiau fc Turck-Chiezei 
2002) the standard stellar evolution model results in too 
high Li depletion during the PMS stage. That is also found 
in my calculations. It is shown in Fig. 3 that the ZAMS 
Li abundance profile of Pleiades predicted by standard 
stellar model(the right dashed line) is too low to fit the 
observations. Applyin g the convection crite r ion modifie d 
by the magne t ic fieldifGough fc Taylerl 1 19661 : iMossI Il968l ). 
IVentura et al.l l|l998bl ) have found that the depletion of 
Li is lower than the value predicted by standard model 
and the difference is determined by the magnetic field in- 
tensity. On the other hand, the stellar activity has been 
thought to lead to the scatter of Li abundance in young 
clust er s JXiong fc Dengll2005l: king et al.ll2010l : iPinsonneaultl 
l20T0h . IXiong fc Dend (|2005l ) have investigated the a Per 
cluster and have pointed out that the dispersion of Li abun- 
dance in the cluster can b e caused by i nhom ogeneous red- 
dening and stellar activity. iKing et alj (|2010l ) have investi- 
gated the cool Pleiades dwarfs and have suggested that the 
striking Li dispersion is caused b y surface activity infl ating 
the star radii in the PMS stage. IPinsonneaultl ((2010) have 
showed that the Pleiades Li dispersion is reproduced when 
the radius of stellar models are inflated 10%. 

The effects of the stellar activity and the magnetic field 
on stellar structure are complex. The complexity makes it 
hard to trace the depletion of Li accurately during the PMS 
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stage. This paper focuses on how the overshooting mixing 
affects Li depletion during MS stage. In order to avoid the 
problem of the PMS Li depletion affected by the stellar ac- 
tivity and the magnetic field, I use the relation of Li abun- 
dance vs. T e ff of the young cluster Pleiades to be the initial 
Li abundance function at ZAMS for all clusters considered in 
this paper. The relation of Li abundance vs. T e / / of Pleiades 
is shown in Fig.3. As mentioned above, the s catter of Li 
abundance could be due to the surface activitv( |King et"al] 
l2010l : |Pinsonneaul3l201Ch . Therefore, the lower envelope of 
Li abundances in Pleiades is set as the initial Li abundance 
of ZAMS stars. The solid line in Fig.3 is the function: 



M/M 



A(Li) = \ S- 1 ' 3 ^ 



;T eff < 6500if 
3.1;T e// > 6500K 



(8) 



which is an approximate profile of the lower envelope of Li 
abundances. Equation (8) is used to calculate the Li abun- 
dance of ZAMS stars in this paper. 



4 NUMERICAL RESULTS 

Using the stellar modeling method described in Section 2, I 
calculated stellar models with different mass in the case of 
six clusters(Hyades, Praesepe, NGC6633, NGC752, M67 & 
NGC3680). Based on those stellar models, I calculated Li 
abundance in both cases with and without the overshooting 
by solving Eq.(2) with the initial Li abundance of ZAMS 
stars calculated via Eq.(8). 



4.1 Theoretical Li abundance isochrones with and 
without the overshooting mixing 

The theoretical Li abundance isochrones of three interme- 
diate clusters(Hyades, Praesepe & NGC6633) are shown in 
Fig.4. The dotted-lines are the initial Li abundance accord- 
ing to Eq.(8). The dashed-lines are the theoretical Li abun- 
dance isochrones of stellar models without overshooting. The 
solid-lines are the theoretical Li abundance isochrones of 
the stellar models with overshooting. It can be found that 
there is almost no Li depletion in those intermediate clus- 
ters in the concerned mass rang when the overshooting is 
absent. However, the overshooting case show modification 
on Li depletion. For lower T E ff star, overshooting mixing 
results in more modification on Li depletion. Lower T e // 
represents higher temperature at the BCE, thus the over- 
shooting mixing is more efficient on depleting Li. The over- 
shooting mixing does deplete more Li, but the modifica- 
tions seems to be not enough to explain the observations of 
clusters Hyades & Praesepe. This is probably for the rea- 
son that the metallicity of Hyades and Praesepe is larger 
than Pleiades's and then Eq.(8) may overestimate the ZAMS 
Li abundance of Hyades and Praesepe. We can estimate 
how much Eq.(8) overestimates the ZAMS Li abundance of 
Hyades and Praesepe via Fig.3. It is found that the predicted 
ZAMS Li abundance of Pleiades is almost identical to Eq.(8) 
at the high temperature rang with T e // > 5500A". This in- 
dicates that the standard stellar model predicts ZAMS Li 
abundance well if T ef! > 5500^. At T e// = 5500A", the 
ZAMS Li abundance profiles predicted by standard model 
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Figure 4. Li abundance of the intermediate clusters Hyades, 
Praesepe & NGC6633. Observation data : circles and triangles 
are the observational Li abundance and the upper limits respec- 
tively. The dotted lines show the relation between Li abundance 
and T^ff at ZAMS. The dashed lines show the Li abundance of 
stellar model without overshooting at corresponding [Fe/H] and 
t, and the solid lines are with overshooting. The upper horizonal 
axis is the mass of stellar model, links to the dashed line and the 
solid one. The upper horizonal axis doesn't link to the dotted line. 
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show about 1.2, 0.9 and 0.8 dex depletion for Hyades, Prae- 
sepe and Pleiades, respectively. Accordingly, Eq.(8) overes- 
timates about 0.4 dex of Li depletion for Hyades and 0.1 
dex for Praesepe at T e // = 5500.ff. If these modifications 
are taken into account in the Hyades and Praesepe cases in 
Fig.4 at T e ff — 55007^, the solid lines should be in better 
agreement with the observations. The gap of Li abundance 
in the rang 6200K < Teff < 6800K in the observations 
of the clusters (Hyades & Praesepe) can't be reproduced by 
the theoretical models because of the absence of the possi- 
ble mechanisms(e.g. atomic diffusion and radiative acceler- 
ation) in this work. It seems to be no Li gap in the cluster 
NGC6633, probably because the observation data in the cor- 
responding rang of T e / / are not enough to reveal the gap. 

The theoretical Li abundance isochrones of three old 
clusters(NGC752, NGC3680 & M67) are shown in Fig.5. 
The observations of Li abundance of NGC3680 and M67 
include only the MS stars with the luminosity lower than 
the turn-off luminosity, in order to distinguish the evolved 
stars which may be in the same rang of T e f / as the investi- 
gated stars. Li abundance of the sun is also denoted in the 
M67 case in Fig.5, since the age of the sun is close to M67's. 
The theoretical Li depletion of the stellar models without 
overshooting are obviously not enough since almost all the 
observation data locate below the theoretical isochrone. It 
is found that the overshooting mixing leads to significant 
modifications on Li abundance in the stellar models. This 
is because the age of those clusters are long, the overshoot- 
ing mixing could deplete more Li comparing with the re- 
sults of the intermediate clusters. The M67 cluster show a 
Li dispersion. It is probably because of the uncertainty of 
T e ff. In deriving the effective temperature of the samples 
of M67 in Jones et al.'s(199 9) wor k , the uniform redden- 
ing value has been adopted. I Taylor] |2007l ) has studied the 
reddening of M67 cluster, and has found that the standard 
error of mean reddening value is a[E(B — V)] = 0.004ma<; 
by analyzing a sample including about 50 stars. Accord- 
ingly, the standard error of the reddening value for a star 
in M67 is a[E(B - V)] = 0.03mag(noting the statistical 
relation a(x) = yNa(x)). For the M67 stars investigated 
in this paper, 0.01 mag uncertainty of color corresponds to 
40K uncertainty of T e ff (Jones et al., 1999). Therefore, the 
standard error of the reddening value for a M67 star leads 
to a(T e ff) = 120K. It can be found in the M67 case in 
Fig.5 that the observations are in agreement with the Li 
vs. T e f f relation predicted by the overshooting stellar mod- 
els, when a(T e ff) — 120K is taken into account. There are 
some stars('I2-a', 'I2-b' and 'S982' in Fig.5) being outside 
la range of Teff, but they are binary components. 

In this paper, the settling of Li is not included because 
the main purpose is to investigate how the overshooting 
affects the Li abundance of the G and K type low mass 
stars. The timescale of the settling is sensitive to the mass 
of star. It is very short for the massive star and is very 
long for the low-mass star. Therefore, the settling should 
not si gnificantly affect the low mass stars. IXiong fc Dend 
(2009) have found that the settling of Li is significant for 
M > I.IMq and is almost no effects for M < 1M . Accord- 
ingly, if the settling is taken into account, the results of the 
low-mass star with T e ff < 5800if in this paper should not 
be changed significantly and the Li abundance of the stars 
with T e ff > 5800if should be depleted more. 
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Figure 5. Similar to Fig.4, but for old clusters NGC752, 
NGC3680 & M67. The solid symbols in the M67 case denote 
binary component. 



4.2 Li burning history and its depletion time scale 

As it is found in Figs. (4-5), the theoretical Li abundance 
isochrones of the stellar models with M ^ 0.90M© show 
only a little Li depletion during the MS stage when the over- 
shooting is absent. In order to investigate the evolution of 
Li abundance of stellar models, the relation between A(Li) 
and t of stellar models with M = 0.90, 1.00, I.IOMq for the 
old cluster M67 are shown in Fig. 6. The solid lines show 
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3.0 




Figure 6. The evolution history of Li abundance of stellar models 
with M = 0.90M Q , l.OOM , I.IOMq for the M67 case. Dash lines 
are for the stellar models without the overshooting and the solid 
lines are for the stellar models with the overshooting. 



the Li abundance of stellar models with overshooting mix- 
ing, and the dashed lines correspond to the no overshooting 
case. For the overshooting case, it is found that the relation 
between A(Li) and t is almost linear. This indicates that 
the surface Li abundance satisfies the equation below in the 
overshooting case: 



[Li] » [Li] exp(- 



t-t , 



(9) 



where tla is the time scale of Li depletion. 

Now, I estimate TLi for the solar mass star. The time 
scale of Li depletion at radius r in the star should be tu (r) = 
Max(TN(r),rov{r)) where rjv(r) is Li burning time scale 
at r and tov(?) is the overshooting mixing time scale. For 
solar mass star, IgTpcE ~ 6.35 during the MS stage, and 
the typical temperature of the reaction Li 7 (p,a)He 4 at the 
location where rjv(r) ~ lGyr is about IgT ~ 6.4(denoting 
this location as A). Therefore, the mixing distance from the 
BCE to A is about L w Q.\2H T ~ 0.3H P (noting that the 
temperature gradient V = dlnT/dlnP m V a d ~ 0.4 near 
the BCE) where Ht = —dr/dlnT. The time scale of the 
overshooting mixing process at A is tov(A) ~ L 2 /Dov ~ 
0.3 2 H P /(CxVk) where Hp ~ 10 10 cm and Vk ~ 10 3 cm/s 
near the BCE in the solar case, thus tov(A) ~ 10 16 s ~ 
lGyr and tn{A) ~ lGyr. tov(t) > lGyr in the region 
IgT > 6.4 and rjv(f) > lGyr in the region IgT < 6.4 because 
the reaction rate R oc T 22 . Therefore, the time scale of Li 
depletion of a solar mass star in the MS stage should be 
TLi = Min{r Li {r)) = r Li (A) ~ lGyr. 

It can be found in Fig. 6 that [Li]/[L£]o ~ 1/e at t fa 
1.5Gyr for the LOOM© star. This means TLi ~ 1.5Gyr for 
the LOOM© star and validates the estimate above. 



5 CONCLUSIONS AND DISCUSSIONS 

I have investigated Li abundance in six open clus- 
ters(Hyades, Praesepe, NGC6633, NGC752, NGC3680 & 
M67) by using the diffusive overshooting approach, which 
describes the overshooting mixing as a diffusion process. 
The diffusion velocity is based on the turbulent convection 



model (|Li fc Yand2007l ; lzhang fc Lil2012bT ). The PMS Li de- 
pletion is affected by the stellar activity and the magnetic 
field. In order to take into account the PMS Li depletion 
and avoid the difficulty of stellar activity and the magnetic 
field, the relation between Li abundance and T e ff of the 
young cluster Pleiades is set as the initial function of Li 
abundance vs. T e // for ZAMS stars. It is found that the 
diffusive overshooting mixing, which is based on the TCM 
with parameters favored by the sun, leads to remarkable Li 
depletion in MS stage for the old clusters. 

There is almost no Li depletion for the stars with 
M ^ 0.9A/q during the MS stage when the overshooting 
is absent. The time scale of Li depletion due to the over- 
shooting mixing is about TLi ~ lGyr for the IMq star in 
the MS stage. For the young clusters(i < lGyr, e.g., Hyades, 
Praesepe, NGC6633), the modifications of the overshooting 
mixing on Li abundance are slight due to t < TLi- How- 
ever, for the old clusters(t > lGyr, e.g., NGC752, M67 & 
NGC3680), the overshooting mixing should be included in 
order to fit the observations. 

A main distinction between this work and Xiong & 
Deng's(2009) work is whether the PMS depletion of Li is 
taken into account. Xiong & Deng's(2009) results show that 
Li is depleted to A(Li) = at T eff « 5200^ for Hyades 
and Praesepe, at T eff w 5500i\" for NGC752, and at T eff « 
5700^ for M67. Comparing Figs. (4) and (5) with Xiong & 
Deng's(2009) results, which ignore the PMS Li depletion, the 
results in this paper show less Li depletion. These indicate 
that the overshooting mixing effects in this paper is much 
weaker than theirs. This is probably because of the small 
parameter(e.g., diffusion parameter Cx) in this model. The 
small parameter is due to the unknown characteristic scale 
of the overshooting mixing process, which should be in the 
range between the Kolmogorov scale and the largest scale in 
the convection zone, and the assumption that the time scale 
of the overshooting mixing sh ould be compara ble with the 
stellar ev olutionary time sca le (|Deng et al]|l996t ). Our previ- 
ous work (|Zhang fc Lill2012ah showed that this weak mixing 
seems to be in consistent with the helioseismic data. 

Noting that the overshooting mixing with the param- 
eters favored by the sun is also suitable for Li abun- 
dance of low mass stars in old clusters, it is thought to 
be reaso nable that the para meters of the TCM in the so- 
lar case l|Zhang fc L3 l2012aTl and the diffusion pa rameter 
Cx ~ 10~ 10 (|Deng et al.lll996l ; IZhang fc Lill2012al ) is gen- 
erally applicable on the downward overshooting of the con- 
vective envelop of low mass stars. 

On the Li abundance in clusters, the overshooting mix- 
ing is not the only candidate. There are other mechanisms 
can deplete Li as mentioned in Section 1. However, t he he- 
lioseismic research l|Christensen-Dalsgaard et al.|[201lT ) have 
showed that the overshooting of the heat flux could pene- 
trate 0.37H P to the location where T = 2.5 x 10 6 K in the 
solar case. And a c cording to s imulations llSingh et al. 19951; 
i Meakin fc Arnettl 120071. l201Ch and TCMs l|Xiong fc Dend 
l200ll ; IZhang fc Lj|2012bh . the overshooting of turbulent ki- 
netic energy, which represents the turbulent r.m.s. velocity 
and affects the mixing directly, can penetrate deeper than 
heat flux. This indicates that the overshooting mixing seems 
to be unavoidable in investigating Li abundance. 
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